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1 Introduction 

T-HIn this paper I will present a new model of eye control 
O developed within the Probabilistic Population Code (PPC) 
framework. The posterior parietal cortex is believed to di- 
^rect eye movements, especially in regards to target track- 
^^ing tasks, and a number of debates exist over the precise 
nature of the computations performed by the parietal cor- 
l ^ tex, with each side supported by different sets of biological 
evidence. In this paper I will present my model which navi- 

i gates a course between some of these debates, towards the 

fljend of presenting a model which can explain some of the 
competing interpretations among the data sets. In particu- 
^lar, rather than assuming that proprioception or efference 
O copies form the key source of information for computing 
eye position information, I use a biological plausible imple- 
mentation of a Kalman filter to optimally combine the two 
signals, and a simple gain control mechanism in order to 
accommodate the latency of the proprioceptive signal. Fit- 
^"ting within the Bayesian brain hypothesis, the result is a 
Bayes optimal solution to the eye control problem, with a 
• range of data supporting claims of biological plausibility. 

1 2 Computation in the Parietal Cortex 

> 

Contemporary research has shown that the posterior pari- 
^ etal cortex is involved in target reaching. In order to accom- 
C3 plish this, the parietal cortex must integrate multi sensory 
information to derive target and body positions, coordinate 
eye movements to track targets, and finally coordinate mo- 
tor activity to reach toward targets. 

The computational characteristics of the parietal cor- 
tex are well characterized and the subject of ongoing re- 
search. Neurons in the parietal cortex exhibit spatial re- 
ceptive fields selective for target and body positions in var- 
ious coordinate systems. Retinal coordinates have been ob- 
served to be especially common in the parietal cortex, and 
eye position information is thought to be required to trans- 
form variables of interest into and out of a standard retinal 
coordinate system (Cohen et al. 2002). The source and 
nature of this eye position information is subject to debate, 
with some arguing for static eye position and others for dif- 



ferential eye position information respectively as the key 
factor in performing these transformations ( jBlohm et al.[ 
2"009t P^uetaLj [20TT] |Wang and Pan| [2012] ) . 



In general, population activity can be partially character- 
ized by its gain - a scaling up or down of the total output 
of the network. Gain normally indicates the precision of an 
encoding in a population, and maintaining gain at various 
stages of neural processing involves respecting the preci- 
sion information of upstream populations. The nonlinear 
operation known as divisive normalization - which gener- 
ally involves dividing a product of gains by their sum - has 
been proven to possess certain powerful optimality proper- 
ties, and has been found to be frequently applied in the low 
level sensory processes of the brain ( |Schwartz et al.[ 2001 1. 

Parietal neurons also exhibit a phenomenon known as 
gain field modulation in which gain encodes some other 
variable of interest such as eye position rather than the 



usual precision information (Blohm and Crawford 2009). 
This mechanism of gain field modulation may also under- 
lie powerful and efficient means of computing coordinate 
transformations (Pouget and Snyder[ |2000[ Chang et al. 



2009), indicating that more than simply being a novel way 
of encoding information, gain field modulation is a critical 
aspect of certain neural computations. 



3 Probabilistic Population Codes 

Within the probabilistic population code framework, a pop- 
ulation of neurons is interpreted as encoding a posterior 
distribution over a stimulus through the stochastic response 
of that population to the stimulus. In this part of the re- 
view I will introduce the basic case of linear, 1 -dimensional 
PPCs. This technical review has been drawn from the arti- 



cles|Pouget an d Snyder] ( |2000p ; |Ma et al.| ( |2006| ) ; |Beck et al. 
(2011 1, and the simulations and plots were developed by 
myself. This review is intended only to provide sufficient 
detail to follow the model presented later. Readers wishing 
to simulate these models themselves are advised to consult 
the aforementioned papers. 

In general, computation in a PPC involves sampling a 
population's response r to a stimulus s, the distribution 
of which is assumed to follow some likelihood distribution 
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p(r | s). A Linear PPC is one in which this likelihood dis- 
tribution is limited to be one from the exponential family 
with linear sufficient statistics - that is, an exponential fam- 
ily distribution where the kernel h is a linear function of s 
0. This equation for the likelihood also contains a multi- 
plier 4>{y) which corresponds to the gain on the encoding, 
through which we can incorporate precision information or 
gain field modulation. 



p{v 



(r)e 



h(s)*r 



(l) 



Given an r sampled from this likelihood distribution and 
a prior over the stimulus a(s), we may decode a posterior 
distribution p(s | r) over s up to a normalizing term Z 
Q. With this distribution we can estimate a value for the 
stimulus encoded by a particular population's response. 



p(s | r) 



h(s)*r+a(s) 

z 



(2) 



Although research has demonstrated that neural activ- 
ity is not always well modelled by Poisson processes with 
stimulus dependant rates, Poisson models remain useful 
and essential as first approximation likelihood models. For 
the simplicity of the mathematics, Gaussian posteriors also 
tend to be easier to work with and are sufficient for many 
purposes. If we assume Gaussian tuning curves around 
some preferred stimulus as defining the rate of our Poisson 
likelihood distributions p(r \ s), then the posterior distri- 
bution p(s | r) takes a Gaussian form. For the rest of this 
review, we will assume Poisson distributed likelihoods and 
Gaussian distributed posteriors. For the case of the static 
encoding of stimuli, this basic formulation is illustrated in 
figure [I] 

3.1 Neural Circuits 

Deterministic computation typically involves algorithms 
which compute the output of a function y = /(x). More 
generally, however, it may be that the dependence of an 
output y on an input x is probabilistic rather than deter- 
ministic, and that what we therefore wish to compute as 
output is a posterior distribution p(y | x) rather than a sin- 
gle value. More general still, instead of being given true 
values of the input x we might be given only distributions 
p(x) over these true values. In this case we must marginal- 
ize over these input distributions in order to calculate an 
output distribution p(y) This marginalization will be 
referred to as a probabilistic computation. 



p(y) 



p(y I x)p(x)dx 



(3) 



The PPC framework is general enough to model neural 
circuits - neural populations which implement probabilis- 
tic computations. Neural implementations of probabilistic 
computations, however, are stochastic and do not explicitly 
encode probability distributions. This means that given a 
PPC and a sample encoding r in of some distribution over 
x, the probabilistic computation implemented by the given 
PPC is computed implicitly over the given encoding to pro- 
duce a sample encoding r out of this posterior distribution 



Poisson PPC with Gaussian Tuning Curves 

O 




(Preferred) Stimulus 
— Tuning Curves OOOMeuron Responses 

Decoded Posterior Distribution; p(s|r) 



Stimulus Posterior Distribution 

Figure 1: Static Poisson PPC Response: Top) A sample response of a 
PPC from the likelihood p(r | s) given the stimulus. Individual neuron 
responses are arranged over the x axis in accordance with their respective 
preferred stimuli. The Gaussian tuning curves are also illustrated. Bot- 
tom) The decoded Gaussian posterior distribution over the stimulus given 
the response from the above plot. 



p(y | x). Moreover, since encodings are sampled from like- 
lihood distributions, there is variability in the encoding of 
the posterior distributions over y. 

In general, solving marginalization equations involves 
solving high dimensional integrals, which is rarely a trivial 
task. Nevertheless, within the PPC framework Bayes op- 
timal likelihood distributions p(r out \ r in ) over encodings 
have been derived for certain interesting cases, providing a 
basic set of neural circuits which can be composed to model 
more powerful computations. 

The gain of a neural population is a simple multiplier on 
population activity, and a key aspect of marginalization in 
neural circuits involves maintaining appropriate gain in the 
output encoding. Since larger responses improve the signal 
to noise ratio, increasing the gain has the effect of increas- 
ing the precision of the encoded posterior. In parallel with 
the biological research mentioned in section [2] it can be 
shown that certain optimal computations within the PPC 
framework implicitly perform divisive normalization, thus 
supporting the biological plausibility of the PPC approach. 

Nevertheless, as also indicated in section[2] gain need not 
encode precision alone, and so optimal precision weight- 
ing need not be the goal in computing a response to input 
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activities. In this case other variables of interest may be 
encoded in the sampled responses from equation [T] via the 
term 4>(r), and downstream populations may later decode 
and use the gain modulating signal as required. Deciding 
what the gain of a PPC encodes is a critical part of apply- 
ing PPCs, and requires careful consideration of the system 
being modelled in order to be made correctly. 



3.2 Linear Coordinate Transformations 

Suppose that we're given the position in head centred co- 
ordinates of a target xr and 'the eye' e# encoded as ap- 
propriate PPC likelihood samples x# and e#, and we wish 
to encode the position of the target in eye centred (reti- 
nal) coordinates x R as a new response xjj. We must then 
solve the probabilistic computation p{A R \ Sc H ,e H ) where 
x R — xh+eh- In order to optimize the performance of the 
computation, the computation must minimize the error in 
an optimality equation over the posteriors Q. 



p(xr I =p(x R | ± H ,e H ) 



(4) 



This Bayesian optimality constraint equates the amount 
of information we have about x R before and after the com- 
putation, thus ensuring that a minimum of information is 
lost in performing the computation. The results of simulat- 
ing a computation satisfying the aforementioned conditions 
are illustrated in figure [2] The gain g 3 of the optimal output 
encoding expressed as a function of the input gains g\ and 
#2 is given in equation [5] This is the expression for divisive 
normalization with two inputs. 
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3.3 Kalman Filters 

Kalman filters are a class of powerful algorithms for com- 
bining observations and model predictions in an optimal 
way. A 1 -dimensional linear Kalman filter is the simplest 
form of such a filter, but is still capable of solving interest- 
ing problems. Where x represents the state and a the model 
rate, u represents the control signal and b the control rate, 
and rj represents the noise, the equation for underlying dy- 
namics of such a Kalman filter has the form of Q. 



dx(t) 
dt 



a * x(t) + b * u(t) + r](t) 



(6) 



Given the parameters of this equation a, b, and rj, a con- 
trol signal u(t), and an estimate of the current state x(t), 
the Kalman filter will make an a priori estimate of x(t + At) 
based on the model dynamics. Then, given a noisy observa- 
tion z(t + At), the Kalman filter will combine this with the 
a priori estimate to produce an a posteriori estimate which 
optimally combines the two signals weighted by their vari- 
ances. Where a is equal to some diffusion parameter A 
and b is 0, equation ([6]) models a simple diffusion process. 
Given this model and a series of observations, a Kalman 
filter can infer the time evolution of this process. 
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Figure 2: Linear Transformation PPC: Top) Sample responses from two 
different populations responding to two distinct stimuli. The summed 
stimulus is indicated by the darker red. Bottom) The decoded posteriors 
from the input populations, and the decoded posterior of the linear co- 
ordinate transformation PPC which attempts to calculate the sum of the 
two input populations. The gain (height) of the summing PPC posterior is 
based on those of the input PPCs via divisive normalization. 



The algorithm for calculating a Kalman filter update step 
involves solving a series of equations, and the posterior dis- 
tribution of the probabilistic computation which respects 
these Kalman equations is given by p(it(t + At) | x(<), z(t + 
At)). It is possible for a PPC to implement this probabilistic 
computation, satisfying an optimality constraint analogous 
to that given in equation [4] If we interpret the rates of 
the Poisson likelihoods from the beginning of this section 
as rates for a set of Poisson processes, we can model a PPC 
over time as a set of stochastic differential equations which 
define the firing rates of the individual neurons of the pop- 
ulation @. 



~dt 



W • x + U ■ x + M ■ 



x • Q ■ x 



(7) 



In this equation, the encoding x of the posterior distribu- 
tion over the stimulus p(s | x) evolves over time. In accor- 
dance with equation [6] W implements the model rate, U 
implements the control rate, and M implements the contri- 
bution of observation. 

The purpose of Q is to act as a regularizer on the time 
evolution of the firing rates, suppressing firing rates when 
they are too high. The Q term in equation[7]is negative and 
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Figure 3: Kalman PPC Simulation: Top Left) The input population gain. In a vision task, low gain can be interpreted as low contrast which causes 
the object to become difficult to see. Bottom Left) Decoding the input population over time. When the gain is low, the input population response 
cannot reliably encode the input. This model is equivalent to the model in (figure [TJ replicated over time. Top Right) Spiking statistics of the Kalman 
population. Bottom Right) Decoding the Kalman population over time. The Kalman population provides a good estimate of the input signal and can 
interpolate values when the input population isn't active. 



quadratic in the population activity x. In simulating the set 
of differential equations [7] the Q term scales the activity 
in proportion to a quadratic function of the activity, thus 
implementing a form of divisive normalization. The results 
from my simulation of a Kalman PPC on a simple diffusion 
process are given in figure [3] 

A few interesting observations can be made of the results 
of this simulation. In figure [5] the sinusoidal gain depicted 
in the top left plot results in distinct yet related effects in 
each other plot. When the gain is low, the input popula- 
tion fires too rarely to reliably encode the stimulus, and 
attempting to decode a posterior in these cases results in 
numerical errors (Bottom Left). Nevertheless, even when 
few spikes are being transmitted from the input population 
to the Kalman population, the Kalman population is able 
to interpolate the missing values (Bottom Right). At the 
same time, a lack of observations causes the standard de- 
viation of the posterior distribution to increase, indicating 
a declining confidence in the estimate of the state of the 
system. Interestingly, the population response statistics of 
the Kalman population show that regions of low confidence 
correspond to higher firing rates (Top Right). However, 
as the firing rates of the individual neurons decrease, the 
range of firing rates exhibited by the population increases. 



Thus, the Kalman population improves its confidence not 
by increasing its mean firing rate, but by maximizing the 
differences amongst the activities of its population. 



The top right plot also indicates that the neurons in the 
Kalman population exhibit exceptionally high firing rates. 
These firing rates are not biologically realistic, but are re- 
quired in order to achieve a high enough signal to noise 
ratio. This could be more realistically modelled by simulat- 
ing a much larger population of neurons, but simulating a 
realistic number of neurons is difficult on a home computer 
since the calculations depend on matrices which scale in 
size to the square of the number of neurons. Nevertheless, 
it's not unreasonable to interpret each neuron in the sim- 
ulation as representing an ensemble of neurons, thereby 
supporting a more realistic interpretation. 

Another problem with the model as implemented is its 
relatively inefficient use of spike information. Spike infor- 
mation efficiency could be optimized based on other com- 
putational approaches ( |Deneve| 2008[ Hennequin et al.[ 
2010p, providing another way to improve the signal to 



noise ratio and allowing us to lower the firing rate of the 
PPC. 
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4 Recent Empirical Findings 

The goal of this paper is to present a biologically realis- 
tic model of eye tracking of targets within the PPC frame- 
work. The structure of the model has been developed by 
surveying the current literature and attempting to navigate 
a path of least controversy through contemporary debates. 
The primary issues of debate within the literature concern 
the relative strengths of eye position signals deriving from 
proprioception versus those deriving from efference copies, 
and role of gain field modulation in neural computation. In 
this section I will review some results from these debates 
which have shaped my model. 

4. 1 Efference Copies vs Proprioception 

In order to implement coordinate transformations, the pari- 
etal cortex requires information about eye position. The 
two primary candidates for being sources of this informa- 
tion are proprioception and efference copies. Propriocep- 
tion is the sense of the relative positions of various parts of 
the body, and the signals it produces derive from various 
sources which monitor the states of muscles, tendons, and 
skin. Efference copies are corollary discharges produced 
by the transmission of motor signals. Efference copies are 
usually construed as copies of the motor commands them- 
selves, and are only transformed into signals pertaining to 
body positions via a forward model of the effect of the mo- 
tor command. 

Proprioception has the capacity to produce static (as op- 
posed to differential) information about eye position. Nev- 
ertheless, it is often assumed to play a role only in the fine 
tuning of an efference copy based system. The primary rea- 
son for this is the apparent delay from the time of motor ac- 
tion to the proprioceptive feedback. For signals measured 
in the somatosensory cortex, contemporary research indi- 



cates the delay to be between 60 and 100 ms (Wang eFaTj 
2007; pCu et al.| 20111. Nevertheless, saccades are usually 
made only about 3 times per second, and between each sac- 
cade the eye tends to remain still. Thus, somewhere from 
60 to 100 ms after a saccade, the delayed proprioceptive 
information remains still and becomes shift invariant, and 
can therefore be treated as accurate and contemporary 

Since motor commands are relative to the current state 
of the system being controlled, efference copies provide a 
differential representation of state. After being fed through 
the forward model, this representation can be integrated 
to form a static representation of eye position. To the ex- 
tent that this signal is fast, reliable, and the forward model 
can be well tuned, it is rarely argued that efference copies 
do not play some important role in providing eye position 
information. The primary debate rather revolves around 
whether the proprioceptive signal is combined online with 
the efference copy to form an estimate of eye position, or 
whether the efference copy is used to derive the online es- 
timate, with proprioception playing only a post estimate 
feedback role ( Ziesche and HamkeiJ |2011[ |Balslev et al. 
20T2||Waiig and Pan||2012]). 



In spite of some desirable properties, the differential, for- 
ward model based position estimate provided by the effer- 



ence copy system suffers from two important draw backs. 
Firstly, some argue that rather than there needing to be a 
distinct integrator of the differential efference copy repre- 
sentation, the differential information could be integrated 
directly into variable encodings themselves. The problem 
with this approach is that in dealing with discontinuous 
stimuli such as blinking targets, when the representation 
of a stimulus temporarily breaks down, the constant of in- 
tegration of the eye position is lost as well. Thus, after 
a blink the differential information is no longer sufficient 
until a resetting of the constant of integration takes place. 
Secondly, even where we independently integrate the dif- 
ferential signal, without feedback the standard deviation 
of a purely forward model based approach eventually di- 
verges. 



4.2 Precision vs Gain Field Modulation 

As described earlier, we have two uses for the gain of a 
PPC - we can either use it to encode precision information 
about the encoded posterior, or use it to encode informa- 
tion about a second variable of interest. In the first case, 
we rely on divisive normalization in order to optimize the 
precision of the encodings of downstream populations. In 
the second case, we discard the precision information and 
effectively turn a PPC which computes a Gaussian posterior 
into a 2-dimensional encoder. 

Previously I had characterized the relationship between 
the choice of precision versus gain field modulation as 
something which is subject to debate in the literature, but 
this is not actually the case. Rather, the literatures study- 
ing optimal encoding and gain field modulation are dis- 
tinct, and when developing a neural circuit there is little 
information about the relative merits of these two encod- 
ing regimes for deciding which regime is appropriate. 

As jChang et~aT] ( |2009P and |Blohm| ( |2012| ) show, there are 
different ways to use gain field modulation to accomplish 
complex transformations. A key characteristic which both 
of these studies share is that the targets which they are 
modelling compute transformations the outputs of which 
are given to the motor system. That is, the outputs of the 
studied neural circuits occur where there is the need to 
compute an actual estimate of the stimulus. In the case of 
Gaussian posteriors, if we use something such as the maxi- 
mum a posteriori (MAP) estimate to evaluate our stimulus, 
the estimate of the stimulus will simply be the mean. Thus, 
the variance of the encoded Gaussian is free to be used for 
some other aspect of the probabilistic computation. 

However, when the response of a given neural circuit is 
likely to be synthesized with other responses downstream, 
it is in most cases necessary to retain precision information 
in order to sensibly weigh the inputs and thus compute the 
optimal output of the neural circuit. This dichotomy is re- 
flected in the literature as well, where models implement- 
ing gain field modulation tend to be of circuits in the pari- 
etal cortex, and those implementing divisive normalization 
tend to be of low level sensory systems. 
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5 Model 

In Bayesian statistics, a recursive Bayesian estimator - of 
which the Kalman filter is a simple form - is an algorithm 
for combining a forward estimate with an observation. 
Where the efference copy system drives the control rate 
matrix, and proprioception provides our observations, the 



r 



Eye 



problem delineated in section 4.1 of estimating eye position 



with proprioceptive and efference copies can be formulated 
as a problem of recursive Bayesian estimation, and poten- 



tially solved by the neural circuit defined in section 3.3 



In translating the delineated problem into one which can 
be solved by a Kalman PPC, one still needs to address the 
question of what to do about the delayed proprioceptive 
signal, as a Kalman filter generally assumes that the obser- 
vation corresponds to the same time step as the forward 
prediction. 

A low gain results in a low precision observation, which 
causes a Kalman filter to rely on its forward prediction. 
As explained previously, the delayed proprioceptive signal 
is accurate from at least 100ms after a motor command. 
Thus, if the efference copy system can also be used to in- 
hibit the response the proprioceptive signal for 100ms after 
any non zero motor command, the proprioceptive signal 
will be inhibited when it is delayed, used when it is effec- 
tively contemporary, and thus combined with the forward 
prediction at exactly those times when it is effective to do 
so. 

Omitting the tildes indicating encodings, the neural cir- 
cuit which I have simulated is depicted in figure [4] We 
begin with the movement of the eye caused by the mo- 
tor command u(t). Where d represents the length of the 
proprioception response delay, eye movement causes pro- 
prioception to transmit a delayed posterior over eye posi- 
tion in head centred coordinates p{eu{t — d)) modulated 
by the gain g(u(t)) transmitted to it from the motor proces- 
sor/ efference copy source. The Kalman PPC integrates this 
gain modulated signal with the forward model driven by 
u(t) and the PPC's internal activity. The Kalman PPC then 
passes its estimate of the contemporary posterior p{eH{t)) 
to the motor processor. Given some target information in 
head centred coordinates xh (e.g. audio target informa- 
tion), the motor processor combines xh with the MAP of 
p(eij(t)) to produce a new control signal u(t). 

The dashed output from the Kalman PPC points towards 
other uses of the eye position estimate. The phenomenon 
of gain modulation has so far been observed primarily as 
supporting the computation of complex coordinate trans- 
formations. Since no coordinate transformations are in fact 
computed in my simple model, sophisticated gain modu- 
lated computation of this kind is not required. Neverthe- 
less, eH is an important signal often encoded via gain field 
modulation, e.g. as the variance of a posterior over some 
target position in retinal coordinates xr, and this output 
signal could be used effect this gain modulation. In this 
vein, a more general version of the PPC reach model imple- 
mented in |Beck et al. ( [2011] ) could be driven by the signal 
en- Through this model gain field modulation could be 
explored within the PPC framework, by observing how en 
affects/modulates downstream populations. 
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Figure 4: Eye Control Model: Representation of the eye control model. 
The dashed boxes represent signal sources, the solid boxes represent PPCs, 
and the semi dashed boxes represent systems which are neither simple 
sources nor PPCs. Target and eye position signals are represented by x 
and e, control and gain signals by u and g, and head centred coordinate 
systems are indicated by the subscript and H. 



The simple saccadic eye control problem solved by my 
model evolves in the following way. Every 0.3 seconds 
a new stimulus position is drawn uniformly from the set 
{—2, —1, 0, 1, 2}. The goal of the model is to track the tar- 
get, and so given a target position and an estimate of eye 
position, the motor processor moves the eye at its maxi- 
mum speed to minimize the difference between the target 
and eye positions. In order to suppress oscillatory feedback 
loops, a motor command is not generated when the differ- 
ence calculated is smaller than some e. 

The Kalman filter is a recursive algorithm, and so in or- 
der to ensure good performance, it is important that the 
internal network dynamics of a Kalman PPC be initialized 
with a stable and precise estimate of the contemporary sys- 
tem state. As such, simulations are begun with a 2 second 
initialization period where the stimulus, eye position, and 
control signal are set to in order to allow the Kalman PPC 
to accumulate evidence for its internal model, thus allow- 
ing transients to expire before we apply the Kalman PPC to 
problem in question. 



6 Results 

The results of my simulations of the model presented in the 
previous section are plotted in figure [5] 

The first plot depicts the gain used to scale the activity of 
the PPC encoding the eye proprioception signal. The time 
evolution of the gain begins with the initialization period 
where it remains at its maximum value. After initialization, 
whenever a non zero control signal is sent to the motor pro- 
cessor, the gain of the proprioceptive PPC is set to near 
(a gain of is not possible as Poisson distributions require 
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Figure 5: Eye Control Simulation: 1) The gain used to scale the activity 
of the proprioception PPC. 2) The decoding over time of the propriocep- 
tion PPC resulting from a simulation of the model. 3) The decoding over 
time of the Kalman population. The Kalman population successfully infers 
the position of the eye and drives the solution to the control problem. 4) 
In different case in which after initialization, the proprioceptive signal is 
withheld. The decoding eventually diverges. 



non zero firing rates). This implements the previously ex- 
plained gain field modulation gating mechanism. 

The second plot depicts the simulated proprioception 
PPC as driven by the gain in the first plot and the state 
of the eye. Close inspection reveals that after every jump in 
eye position, a range exists where the eye position cannot 
be decoded from the proprioception PPC. The purpose of 
this gain modulation is to filter out those parts of the pro- 
prioceptive encoding where the delayed and contemporary 
signals differ. Thus, all that remains in the activity of the 
PPC is an encoding of the contemporary position. 

The third plot depicts the successful inference of eye po- 
sition, and the successful tracking of the target stimulus 
with the eye. The initialization period drives the standard 
deviation of the encoding down to a minimum, providing 
the PPC with enough information to make accurate forward 
predictions. After the initialization we have only intermit- 
tent access to the proprioceptive signal. The regular drop- 
ping out of the proprioceptive signal results in little effect 
on the quality of the estimate, and the estimate is accurate 
enough to drive the motor processor to successfully track 
the target stimulus. 

the final 



Corroborating the claim made in section 4.1 



plot of figure [5] displays a modified version of the com- 
plete simulation. In this case we initialize the network as 
per usual, but then set the gain of the proprioception PPC 
nearly to zero for the rest of the simulation. The result, 
as expected, is the slow divergence of the prediction, and 
eventually divergence of eye control itself. By the end of 
the simulation the tracking of the target has almost com- 
pletely decohered. 



7 Discussion 

In summary, in this paper I have presented a PPC neural cir- 
cuit which efficiently solves the saccadic eye control prob- 
lem while satisfying certain mathematical and biological 
constraints. Mathematically, my model optimizes the com- 
bining of proprioceptive signals with efference copy signals 
through a combination of a Kalman filter with a gain modu- 
lation based gating mechanism. Biologically the model fits 
well the observations made in section [4} and theoretically 
the model is implemented as a set of Poisson processes im- 
plementing Bayes optimal computations, setting it within 
established literature on modelling neural activity by Pois- 
son processes, and the Bayesian brain hypothesis. 

Further work pointed at by this paper include the sim- 
ulation of higher dimensional models which rely on non- 
linear coordinate transformations. In the simple model 
which I simulated, no coordinate transformations needed 
to be computed, and thus gain field modulation was not 
required to compute solutions to, e.g. higher dimensional 
non linear coordinate transformations. Nevertheless, it is 
clear that this is due only to the simplicity of the model, 
and that a more complex model would require that we in- 
vestigate computations which would likely necessitate gain 
field modulation. At the same time, groundwork has been 
laid for exploring various kinds of gain modulated com- 
putation using probabilistic population codes, indicating 
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that this question of the computational necessity of gain 
field modulation could be well investigated within the PPC 
framework. 

The use of gain as a gated signal provides another way 
we might look at gain field modulation and the kinds 
of computation it may help implement. Moreover, find- 
ings in |Xu et alT] ( |201ip and unpublished findings from 
Xu, Karachi, and Goldberg indicate that immediately af- 
ter a motor command but prior to the onset of an accu- 
rate proprioceptive signal, proprioceptive representations 
of eye position in the somatosensory cortex are somehow 
disrupted. This preliminary evidence provides support for 
the biological plausibility of the gating mechanism imple- 
mented here, in particular for the proposal that there might 
be an interaction between corollary discharges and propri- 
oception. I hope that further investigations into these pre- 
liminary experimental findings might be assisted by the re- 
sults of this paper, towards the goal of bringing clarity to 
our understanding of neural computation in the parietal 
cortex. 
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